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We calculate next-to-leading-order (NLO) corrections to the B -> i transition form factors at 
leading twist in the kr factorization theorem. Light partons off-shell by k% are considered in the 
quark diagrams, in the effective diagrams for the B meson wave function defined with the effective 
heavy-quark field, and in the effective diagrams for the pion wave function. It is explicitly demon- 
strated that the infrared logarithms In k% cancel between the above sets of diagrams, as deriving 
the fcT-dependent NLO hard kernel from their difference. The infrared finiteness of the hard kernel 
confirms the application of the kr factorization theorem to B meson semileptonic decays. The NLO 
pion wave function is identical to those constructed from the pion transition and electromagnetic 
form factors, consistent with its universality. Choosing the renormalization and factorization scales 
lower than the B meson mass, the NLO corrections are under control: they amount only up to 30% 
of the form factors at large recoil of the pion, when varying models for the meson wave functions. 
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I. INTRODUCTION 



B meson transition form factors are an essential input of the factorization approaches to nonlcptonic two-body 
B meson decays, such as the perturbative QCD (PQCD) approach [J Q based on the kr factorization theorem 0- 
§]. For next-to-leading-order (NLO) contributions in leading-twist PQCD, the vertex corrections, the quark loops, 
and the magnetic penguins associated with the weak decay vertices in factorizablc emission amplitudes have been 
calculated [9l-[il|. As explained in |{|, the above corrections may be the most crucial NLO pieces for understanding 
the known B — ¥ tttt and B — > ttK puzzles, which result from the large observed 7r°7r° branching ratio, and from the 
dramatically different direct CP asymmetries between the ir =F iv" ± and -k^K^ modes, respectively. There have been 
many applications of this NLO PQCD formalism to nonleptonic two-body B and B s meson decays in the literature. 
For NLO corrections to spectator diagrams, we have identified the so-called Glauber divergences, in additional to 
those which are absorbed into hadron wave functions, and summed them into a phase factor to all orders fl2j . It was 
observed that the phase factor modifies the interference pattern between the spectator diagrams, and further improves 
the resolution of the B — > hit, ttK puzzles in NLO PQCD. At the same level of accuracy, we need to calculate NLO 
corrections to the B meson transition form factors for completeness. 

In this paper we shall extend the NLO framework for the pion electromagnetic form factor in the kx factorization 
[l3| to the B — > 7r transition form factors. In this framework light partons in both QCD quark diagrams and effective 
diagrams for hadron wave functions are off mass shell by k\ [lj, [l5| ■ Not only the collinear divergences from gluon 
emissions collimated to the pion, but also the soft divergences from gluon exchanges between the two mesons exist. 
Compared to the pion form factor (T^ |. a new point is that an infrared regulator associated with the b quark is not 
needed. Due to its finite mass, gluons radiated by the b quark do not generate collinear divergences. Soft divergences 
can be regularized cither by the virtuality of internal particles, or by the virtuality k\ of other light partons, to which 
the radiative gluons attach. That is, the b quark remains on-shell in the above framework, a condition which justifies 
the approximation of the b quark field by the effective heavy-quark field for defining the B meson wave function. 
The NLO pion wave function is found to be identical to those constructed in the pion electromagnetic and transition 
form factors [HI, [lj| , consistent with its universality. Note that the diagrams considered here differ from those in the 
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Q CD- improved factorization (QCDF) approach and in the soft-collinear effective theory (SCET) [13], which are 
based on the collinear factorization theorem (l8j : there is no end-point singularity in the kx factorization, so it is not 
necessary to introduce soft form factors [l9| or to perform the zero-bin subtraction [2(| in our calculation. 

It will be demonstrated that the collinear and soft divergences in the quark diagrams are cancelled by those in the 
effective diagrams for the B meson and pion wave functions. Taking the difference of the above sets of diagrams, 
we derive the fcr-dependent NLO hard kernel at leading twist for the B — > n transition form factors. The infrared 
finiteness of the hard kernel confirms the application of the kx factorization theorem to B meson semileptonic decays 
[2lj . Similar to the analysis in (l3l. Il5j. both the large double logarithms a s In 2 kx and a s In 2 x, x being a parton 
momentum fraction, are identified. The former is absorbed into the B meson and pion wave functions and summed 
to all orders in the coupling constant a s by the kx resummation [l[ , and the latter is absorbed into a jet function 
and summed to all orders by the threshold resummation [22]. Due to the dominant soft dynamics associated with 
the b quark, the effect of the kx resummation from the B meson side is minor. The renormalization scale fi and the 
factorization scale /if are introduced by higher-order corrections to the quark diagrams and to the effective diagrams, 
respectively. Choosing /i and /if appropriately, with both being lower than the B meson mass as postulated in [H, 0] , 
the NLO corrections are under control: they amount only up to 30% of the form factors at large recoil of the pion, 
when varying models for the meson wave functions. 

In Sec. II we calculate the 0(a 2 ) QCD quark diagrams for the B — ► -kIv semileptonic decay, the 0(a s ) effective 
diagrams for the B meson and pion wave functions, and their convolutions with the 0(a s ) hard kernel. Since the 
kx factorization is appropriate for QCD processes dominated by contributions from small x [14| |. we shall keep only 
terms in leading power of x. The important double logarithms are identified, and the fc^-dependent NLO hard kernel 
is presented. Section III contains the numerical investigation, in which we examine the dependence of the NLO 
contributions to the B — > it transition form factors on the renormalization and factorization scales, and on the shape 
of the B meson and pion wave functions. Section IV is the conclusion. 



II. NLO CORRECTIONS 



In this section we calculate the 0(a 2 s ) quark diagrams for the B — > -kIv semileptonic decay, and the 0{a s ) effective 
diagrams for the B meson and pion wave functions in the Feynman gauge. The B — 7r transition form factors are 
defined via the matrix element 



(7r(P 2 )|u 7 ^|B(Pi)) = f + tf)(P? + Pit) + I/V) - fV)]' 



(1) 



where raj (fn n ) is the B meson (pion) mass, and q = P± — P2 is the transfer momentum. The momentum Pi 
(P2) of the B meson (pion) is chosen as Pi = Pi + (1, 1,0t) (P2 = (0,P 2 ~,Ot)) with the component P^ = rn B /y/2 
(P 2 ~ = rjm B I 'y/2). The large recoil region of the pion corresponds to the energy fraction rj ~ 0(1). According to the 
kx factorization, the anti-quark q carries the momentum k% = (x\P^~ , 0, kir) in the B meson and k 2 = (0, x 2 P^~ , k 2 r) 
in the pion, x\ and x 2 being the momentum fractions, as labelled in the leading-order (LO) quark diagrams in Fig. [T] 
We postulate the hierarchy 



m 2 B » Xim B > xim 2 B > xix 2 m 2 3 , k\ Tl k 2 T , 



(2) 



in the small- x region, which is roughly consistent with the order of magnitude: x 2 ~ 0.3, x% ~ 0.1, ra B ~ 5 GeV, and 
kx < 1 GeV. Under the above hierarchy, only those terms that do not vanish in the x — > and kx — > limits are 
kept, so the expressions of our NLO results will be greatly simplified. 

To obtain the LO hard kernels, we sandwich Fig. [T] with the following leading- twist spin projectors for the B meson 
and the pion 0, [H 
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is the number of colors. The contributions proportional to the B meson distribution amplitudes and 4>\, ' from 



respectively, where the dimensionless vectors are defined by n + = (1,0, Ox), and n 
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FIG. 1: Leading-order quark diagrams for the B — !> 7r transition form factors with <g> representing the weak vertex. 



with the strong coupling g, and the color factor Cp. To reach the above expression, we have applied the hierarchy 
X2m 2 B 3> k\ T in Eq. @ to the internal b quark propagator. The denominator x\X2"r\vn 2 B + |kiT — k2T| 2 comes from 
the virtuality of the LO hard gluon, in which the |kiT — k2r| 2 term smears the end-point singularity from small X2- 
Similarly, Fig. Hfb) leads to the amplitude 

H y b >{x 1 ,k 1Tl X 2 ,k2T) = C F f— fr- <Mx 2 ). (5) 

r)(xiX2irm B + \ki T - k 2T r) 

Comparing Eqs. (01 and ([5]), it is easy to see that the term proportional to 4> B ' from Fig.QJa) dominates numerically 
according to the hierarchy in Eq. ([2]). As explained above, the B — > n form factors receive major contributions from 
the small- a; region, in which the kj< factorization is an appropriate framework. Since the amplitude from Fig. [TJ(b) is 

suppressed by a power of X2, we will not consider the NLO corrections to H^°\xi, kyr, X2, k^r), and focus on those to 
Fig- da) below. The term proportional to in Eq. ([5]) gives the symmetry breaking effect [lg], which is calculable 

even in the collinear factorization, as convoluted with (j)~ B \xi) ~ x\ at small x\. 



A. NLO Quark Diagrams 

The NLO corrections to Fig. HJa) contain Figs. [3J [31 and 0] for the self-energy corrections, the vertex corrections, 
and the box and pentagon diagrams, respectively. The ultraviolet poles are extracted in the dimensional reduction 
24] in order to avoid the ambiguity from handling the matrix 75. We adopt the following convenient dimensionless 
ratios 
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as presenting our results. The infrared poles arc then identified as the logarithms ln^i and In #2. 
The self-energy corrections in Fig. [2] give 
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where 1/e represents the ultraviolet pole, /i is the renormalization scale, je is the Euler constant, Nf is the number 
of quark flavors, and H^ ' denotes the leading-twist LO hard kernel proportional to P^, 

H^( Xl ,k 1T ,x 2 ,k 2T ) = - 4g2( r FP f . (12) 

x 2 r]di 2 m B 

The above expressions are basically similar to the corresponding ones obtained in the pion electromagnetic form 
factor 13]. We emphasize only that Fig.^a), the self-energy correction to the b quark, requires a mass renormalization 
as indicated by the first term in the square brackets of Eq. ([7]) . The finite piece of the first term is then absorbed, 
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-kf T , into the redefinition the b quark mass, 
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In this work we shall not differentiate mb(p) from tub, because the distinction between them contributes at next-to- 
leading power. The second term in the square brackets of Eq. ([7]) represents the correction to the b quark wave function. 
As explained before, we shall consider an on-shell valence b quark, so the involved soft divergence is regularized by a 
gluon mass m g , which will be cancelled by the corresponding soft divergence in the effective diagram Fig. (Sfa) below. 
The results from the vertex corrections in Fig. [3] are summarized as 
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FIG. 3: Vertex corrections to Fig. [lja). 
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FIG. 4: Box and pentagon diagrams. 



The amplitude from Fig. Eta) depends only on the regulator 82, because the radiative gluon attaches to the virtual b 
quark line. The double logarithm 2 In $2 hi £2 leads to the known Sudakov logarithm In 2 82 and the known threshold 
logarithm In 2 X2 [IH, [l5[ , as reexpressed in the form 



2 In S 2 In x 2 = In 2 5 2 + In 2 x 2 - In 2 — . 

X2 



(20) 



The radiative gluon in Fig.[3tb) attaches to the massive valence b quark and the virtual b quark, so Eq. (I16|) is infrared 
finite. The radiative gluon in Fig. [H^c) attaches to the light valence anti-quarks, such that both the collinear and soft 
divergences are produced, with the latter being denoted by the product ln<5i ln^- This term can be absorbed neither 
into the B meson wave function nor into the pion wave function. Since the radiative gluon attaches to the virtual LO 
hard gluon in Fig. E^d), the soft divergence does not appear Eq. (fT8|) . Equations (fT7|) and (fl"8|) are symmetric under 
the exchange of the regulators Si and 82, as they should. Similar to Fig. I^b), Fig. [3£e) also gives an infrared finite 
contribution. 

The box diagrams and the pentagon diagrams in Fig. |4]lead to the amplitudes 
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Note that Eq. (1211) is power-suppressed in the small x 2 region, while the corresponding diagram gives a leading 
amplitude in the"plon form factor 0. The difference is attributed to the spin projectors: it is A- 

oc 7+ on the B 

meson side here, but f\ oc 7~ on the initial pion side in the latter case. Simply counting the sequence of the gamma 
matrices, it is easy to understand that Fig.SJa) does not produce an amplitude proportional to H^> at leading power. 
Figure HJb) is a two-particle reducible diagram, so its contribution will be cancelled by the corresponding effective 
diagram for the pion wave function [13j , and needs not to be computed. Figure EJc) also contains the soft divergence 
denoted by the Info In fo term, which cancels that in Fig.[3][c). It seems that Fig.EJd) generates a collinear divergence, 
as the gluon on the right is parallel to the light anti-quark in the pion. However, a careful look at the sequence of the 
gamma matrices, similar to that for Fig. Eta), reveals power suppression on this collinear divergence. Equation (|24[) 
does not depend on an infrared regulator associated with the massive valence b quark, because S 2 alone is enough to 
regularize the collinear and soft divergences. The collinear divergence associated with the light valence anti-quark on 
the B meson side is also power-suppressed in Fig. @|f), so Eq. (1251) does not contain Info. 
The amplitudes from all the NLO quark diagrams are summed into 
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for Nf = 6. The ultraviolet pole in the above expression is the same as in the pion electromagnetic form factor, which 
determines the renormalization-group (RG) evolution of the coupling constant a s . 



B. NLO Effective Diagrams 

The 0(a s ) B meson wave function $^ 25] and the 0(a s ) pion wave function [TJ, [H[ collect the effective 
diagrams from the matrix elements of the leading Fock states 
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respectively, with z — (0, z~,%t) an( l V = (y + ;0, yr) being the coordinates of the anti-quark field q, respectively, h v 
the effective heavy-quark field, and T an appropriate gamma matrix. In the above expressions the Wilson line W z {n{) 
with nf is written as 



W z (ni) = Pexp 



ig / dXni ■ A(z + Xni) 



(29) 



and the definition of the Wilson line W v (n 2 ) is similar. It is understood that the two Wilson lines W z (ni) and Wo(ni) 
(W y (n 2 ) and Wo(n 2 )) are connected by a vertical link at infinity jUH^]. Equation (f27|) ( ([28]) ) produces additional 
light-cone singularities (25l. [29l |30]| from the region with a loop momentum collinear to n_ (n+), as the Wilson line 
direction approaches the light cone, i.e., as n\ — > n_ (n 2 — > n + ) (29j . Hence, n\ and n 2 serve as the infrared regulators 
for the light-cone singularities in our formalism. The B meson and pion wave functions then depend on the scales 
£ 2 = 4(ni • Pi) 2 /\n\\ and (f = 4(n 2 • P 2 ) 2 /\n 2 \, respectively, whose variation is regarded as a factorization-scheme 
dependence. This scheme dependence, entering the hard kernel when taking the difference between the quark diagrams 
and the effective diagrams, can be minimized by adhering to fixed n 2 and n 2 . 
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FIG. 5: 0(a s ) diagrams for the B meson wave function. 



We compute the convolution of the NLO wave function & B with the LO hard kernel over the integration 
variables x\ and k'ly, 



$«®£T(°) = j dx[d 2 k\ T ^\x 1 ,k 1 T;x' 1 ,k' 1T )H^(x' 1 ,k' 1T ,x 2 ,k 2 T). 
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The sign of the plus component ri^ of the vector n\ is arbitrary, which could be positive or negative (n 1 has a positive 
sign, the same as of P 2 ~)- Choosing nf < 0, i.e., n\ < as in [H0, H3, we derive, from Figs. |5{a)-|3g), 
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/Xf being the factorization scale. The two-particle reducible diagrams Figs. [5ja) and [Sfc) are calculated, since the 
effective heavy-quark field employed in the B meson wave function differs from the b quark field in the quark diagrams. 
Though the effective diagrams and the quark diagrams have the same soft poles, the finite pieces are different, which 
contribute to the NLO hard kernel. The self-energy corrections to the Wilson lines in Figs. [5jh)-[5jj) yield 
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the same as in the pion wave function 

It is pointed out that the gluon mass m g has been adopted to regularize the soft divergences in the diagrams 
involving the effective heavy-quark field, namely, Figs. Ufa), [5jd), and[5je). The soft divergence in Fig. [5fa) indeed 
cancels that in Fig. [^a) as stated in the previous subsection. The m g dependence disappears in the sum of Eqs. 
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FIG. 6: 0(a s ) diagrams for the pion wave function. 



and (|3"5")l , which must be the case, because the gluons emitted by the b quark and attaching to other particle lines do 
not generate soft divergences. The hierarchy £f >• was employed in the derivation of Eq. (|35|) |25| . so the large 
double logarithm \n 2 (C, 2 / m 2 B ) demands an additional resummation treatment of the B meson wave function, which 
will not be performed in this work. The double logarithms In 2 S± from the quark diagram Fig. E^d) and from the 
effective diagram Fig. Efc) cancel each other. The double logarithms In 2 (m 2 B 8i / (x 2 (, 2 )) are not only attenuated by 
xf, but also cancel exactly between Eqs. (|36| and (|37|) . Summing all the above 0(a s ) contributions, we obtain 
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We then compute the convolution of the NLO wave function $a) with the LO hard kernel over the integration 
variables x' 2 and k'2T, 
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The corrections from Figs. |6][a)-[6]j) are summarized as 
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which are similar to those extracted from the pion transition and electromagnetic form factors [13l . Il5| . but with 
the hard scale Q being replaced by tub here. This similarity supports the universality of the pion wave function. 
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Summing all the above 0(a s ) contributions, we have 
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We stress that the ultraviolet poles are different in Eqs. (|3"9")l and PS)) , since the former involves the effective heavy- 
quark field, instead of the b quark field. That is, the B meson and pion wave functions exhibit different evolutions as 
illustrated below. First, the B meson decay constant, defined via the matrix element with the effective heavy-quark 
field, evolves with an energy scale. Hence, part of the ln/zf term in Eq. (|39j) should be absorbed into /s(/Xf) through 
the RG equation in the heavy quark effective theory (HQET) 
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7/ /s(Mf) = 0, 



(49) 



with the anomalous dimension 7/ = —3 at one loop [33 ] . The RG equation for the B meson wave function without 
the decay constant, (&b(xi, Hf)/ fsiHi), is then written as 



^ dfj, 



a s C F \ $B(xi,[if) 

-IB 



4?r 



0. 



(50) 



where the anomalous dimension 



IB 



1 \ ln^f + 2 



(51) 



governs part of the RG evolution in the kr factorization formulas for the B —> ir form factors p], [2l[ . 



C. NLO Hard Kernel 



The infrared-finite fcr-dependent NLO hard kernel for the B — > tt transition form factors is derived by taking the 
difference between the quark diagrams and the effective diagrams [3] 

-ff (1) (iri,ki T ,:E2,k2T) = G (1) (xi,ki T ,a;2,k2T) 

- J dx'^xT^ {xi , k 1T ; x\ , k' 1T )ff (0) (x[ , k'lT, x 2 , k 2T ) 

-J dx' 2 d 2 k' 2T H (a \x 1 ,k 1T ,x' 2 ,k' 2T )^i 1 \x 2 ,k 2 T;x' 2 ,k'2T). (52) 

Note that a s appearing in Eqs. (|26|). ([39]). and (j48]) denotes the bare coupling constant, which can be rewritten as 

a s = a s (pt) + 8Z{nt)a s {iMi), (53) 

with the counterterm SZ being defined in the modified minimal subtraction scheme. We insert Eq. (|53j) into the 
expressions of the LO and NLO quark diagrams, and of the NLO effective diagrams. The LO hard kernel 
multiplied by SZ then regularizes the ultraviolet pole in Eq. (f2l)|) . The ultraviolet poles in Eqs. (|3"9"|) and (|4*51) are 
regularized by the counterterm of the quark field and by an additive counterterm in the modified minimal subtraction 
scheme. 

The NLO hard kernel for Fig. [lja) is given by 
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in which all the infrared regulators m g , 8%, and <5 2 have disappeared. A choice of the scales Ci and C2 corresponds 
to a factorization scheme, which should be fixed for consistency. Following the scheme ( 2 = Q 2 adopted in the NLO 
analysis of the pion transition and electromagnetic form factors [HI, [l5[ , we set £2 to m 2 B . The important logarithms 
ln(m^/Ci)) arising from the B meson wave function, enter the hard kernel after the infrared subtraction. Instead 
of performing resummation of these logarithms, we choose a sufficiently large (1, sa y, Ci/ m B = 25 in the numerical 
analysis, which has been assumed for achieving the simplified result in Eq. (f35|) . In this scheme the ln 2 (ms/Ci) term 
happens to cancel the large constant term in the hard kernel, and reduces the NLO correction. 

Moreover, the double logarithm In 2 a; 2 has been absorbed into the jet function J(a; 2 ) [22| defined in the kinematic 
region where the virtual b quark in Fig.QJa) becomes almost on-shell, namely, with X2 — > 0. The organization of this 
important logarithm to all orders leads to the threshold resummation factor, which further suppresses the end- poi nt 
singularity from small %i in the B — > n form factors [2lJ. Therefore, we have to subtract the NLO jet function |15| 



J(D#(0) 

from Eq. (|54l) . which finally turns into 



a s (fJ,f)CF 
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(56) 



Because the double logarithm In 2 x\ was not resummed in [2l|, it is left in the above NLO hard kernel Another 
double logarithm In (2:12:2) actually arises from the approximation In 2 <5 12 ~ In (2:12:2). This approximation makes 
sense: a logarithm does not develop an end-point singularity, so the k 2 -, term is negligible in ln<5i2. Equation (|56l) . 
proportional to P2 11 , generates the NLO corrections at leading twist to the B — >• it transition form factors f + (q 2 ) and 
f°{q 2 ) in Eq. ©. 

III. NUMERICAL ANALYSIS 

In this section we evaluate the B — > 7r transition form factors numerically in the kx factorization up to NLO, 
adopting the following non-asymptotic pion distribution amplitudes 'A'A. 34j , 




x(l - x) 1 + a 2 Cl /2 {u) + a 4 C 3 4 /2 (u) 



1 + 0.59C 2 1/2 (u) + 0.09C 4 1/2 (u) 



?=x(l-x) 1 + 0.11C 2 3/2 (m) 



(57) 



2^2N C 

with the pion decay constant f n = 130 MeV, the Gegenbauer moments 02 = 0.16 and 04 = 0.04, and the Gegenbauer 
polynomials 

Cl /2 (u)=U, C 3/2 (u) = 3u, 

3, 



Cl /2 {u) = \{iu 2 - 1), Cl /2 {u) = -(5u 2 1), 



Cl /2 {u) = ±u(5u 2 - 3), Cl'*(u) = ^(21m 4 - 14u 2 + 1) 

8 



15. 



(58) 



and the variable u — 1 — 2x. The B meson distribution amplitudes inspired from a QCD sum rule analysis in the 
HQET [HI 

2 



<4 +) (x,&) 
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0.20 

FIG. 7: Renormalization scale t s (fjb{), defined in Eq. (|61[) . as a function of momentum fractions x\ and X2 for a typical 
factorization scale [if = 1.5 GeV and the ratio Ci/rriB = 25. 



are employed, where the B meson decay constant is set to a constant /s = 214 MeV for convenience (namely, 
neglecting its evolution), and the shape parameter is chosen as loq — 0.35 GeV. 

The first issue concerns the choice of the renormalization scale fi and the factorization scale [if in order to minimize 
the NLO corrections to the form factors. For the first choice, /if is set to the hard scales specified in the PQCD 
approach to exclusive processes @, 0, HH 



f a = ma,x(y/x^rjm B , l/6i, I/&2), t b = ma,x(^/x^rjm B , l/&i, I/&2), 



(60) 



corresponding to the largest energy scales in Figs. Ufa) and[IJb), respectively. Then we utilize the freedom of choosing 
fi to diminish all the single-logarithmic and constant terms in the NLO hard kernel, which is found to be 



t s (fx { ) = \ Exp 



with the coefficients 
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To have an idea of the magnitude of the renormalization scale fi — t s (fi{), we display its behavior in the dominant 
region with the small momentum fractions X\ and x-i in Fig. where the factorization scale [if is fixed at its typical 
value 1.5 GeV. The ratio of the NLO contributions over the total ones as a function of the transfer momentum squared 
q 2 is summarized in Fig. El which is approximately 30% for both the form factors f + (q 2 ) and f°(q 2 )- In the second 
choice, we set both scales to fi = /if = t a (t h ) in the factorization formula associated with Fig.QJa) (JTJb)). It turns out 
that this simple scenario yields larger NLO corrections around 40% as shown in Fig. [SJ The third choice corresponds 
to /if = ms and /i = t s (/if), for which the NLO corrections are approximately 15%. However, the inverse relation 
[if > [i in this case seems not to be natural. Hereafter, we shall adopt the first choice of the renormalization and 
factorization scales as the default one. 

The B — > 7T transition form factors in the kr factorization up to NLO are presented in Fig. [S] It is observed that 
the LO and NLO contributions exhibit the similar power-law behavior, as they should. It is not a surprise that the 
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FIG. 8: Ratios of the NLO corrections over the total contributions to the B — > ty form factors for three different choices of the 
renormalization and factorization scales. 




form factors at the maximal recoil of the pion, / + (0) = /°(0), are close to their LO value [2l|, even after including 
the NLO contributions. The reason is that the meson wave functions have been adjusted accordingly to maintain this 
value, which is regarded as an input. That is, when choosing hadron wave functions in the PQCD approach, one must 
pay attention to the order in the coupling constant, at which the hadron wave functions are determined. Though the 
form- factor values, treated as inputs, are not changed at higher orders, the different hadron wave functions extracted 
at different orders do affect other topologies of nonleptonic two-body B meson decay amplitudes. It is worthwhile to 
investigate the corrections to nonleptonic two-body B meson decays from this NLO source in future works. 

To test the impact of higher conformal-spin partial waves in the pion distribution amplitudes, we plot the q 2 
dependence of the form factors in Fig. [TU] with the asymptotic pion distribution amplitudes. Numerically, both the 
form factors are reduced by about 25% for q 2 < 12 GeV 2 without the non-asymptotic Gegenbauer terms in Eq. (|57l) . 
We also investigate the effects from different models of the B meson distribution amplitudes. A model widely adopted 
in the PQCD analysis is given by 



4 +) M) 



4> ( B \x,b) = 
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■ N, 



■x 2 (l- a:) 2 Exp 



x 2 niB 



(62) 
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FIG. 10: LO and NLO contributions to the B — > ir form factors with the asymptotic pion distribution amplitudes and the first 
scenario for the scale choice, fit = t and n — t s (nf). 




FIG. 11: LO and NLO contributions to the B — >• n form factors with the B meson distribution amplitudes in Eq. (|62[) and the 
first scenario for the scale choice: fj,f — t and fi = t s (jtt/). 



with the normalization constant Nb denned via / dx(j} { g'(x,0) = /s/(2 v / 2]Vc). Note that the two leading B meson 
distribution amplitudes have been assumed to be equal for the purpose of numerical estimate, which do not obey the 
equations of motion (37j . Besides, it exhibits an asymptotic behavior at x — > different from that derived in [38]. The 
corresponding q 2 dependence in Fig. [TT] indicates that the form factors with the model in Eq. (I62|) are approximately 
25% smaller than those with the model in Eq. (|59l) . It is interesting to notice in Fig. [11] that the NLO corrections 
are relatively small, less than 20% of the total contributions. The reason is attributed to the fact that the end-point 
region of x\ is strongly suppressed by this model and the double logarithm In 2 x\ in the NLO hard kernel does not 
play an essential role. 

The extraction of the Cabibbo-Kobayashi-Maskawa matrix element \V u b\ [39( from the semileptonic decay B — > tt£D 
is of intensive phenomenological interest recently (see [34| and references therein). Here we comment on the consistency 
of the B —> 7r form factors predicted in the NLO kx factorization with those in the literature, in view of the extraction 
of \Vub\- For this purpose, we also estimate the theoretical uncertainties of the form factors f + (q 2 ) and f°(q 2 ) from 
the variations of the Gegenbauer moments and 04 in the twist-2 pion distribution amplitudes, from the variations 
of the chiral scale mo = m 2 /(m u + m^), m u (mj) being the u (d) quark mass, involved in the two-parton twist-3 pion 
distribution amplitudes [40( , and from the variations of the shape parameter uiq in the B meson distribution amplitudes 
in Eq. (l59|) . Unfortunately, the extraction of ujq still suffers large uncertainty from QCD sum- rule calculations. We 
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FIG. 12: LO and NLO contributions to the B — ¥ tt form factors with the B meson distribution amplitudes in Eq. (|59[) . however, 
varying the shape parameter loq from 0.30 GeV to 0.40 GeV and the first scenario for the scale choice: fif =t and n — t a (fif). 
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FIG. 13: Theoretical uncertainties of the B — > tv form factors with the first scenario for the scale choice: fif = t and n = t s {fif) 



simply take luq = 0.35 ± 0.05 GeV to illustrate the effect on the form factors from the variation of ojq. It is seen from 
Fig.[TJ]that both the form factors f + (q 2 ) and f°(q 2 ), including LO and NLO contributions, increase (decrease) by 15 % 
with the decrease (increase) of uj . Combining the uncertainties due to 02(lGeV) = 0.16^07' 04(lGeV) = 0.04^q; 08 , 
77i (lGeV) = 1.74±°:3sGeV, and w = 0.35 ± 0.05 GeV, we predict the form factors f+(q 2 ) and f°(q 2 ) as displayed 
in Fig. [13] Fitting to the BaBar data on the integrated B — > tt£D branching ratio within the region < q 2 < 8 GeV 2 
[4lj . where the leading-twist kr factorization is expected to work well, we obtain 

l^|- 2 -90lu'Io| t/ ,t°o:^l e:c ,- (63) 

The above value is in good agreement with that in |4lj . which employed the data on q 2 bins in the whole kinematic 
region and the lattice QCD results of the B — > 7r form factors from the FNAL/MILC Collaboration Equation (|63l) . 
however, differs from \V n },\ = 3.59^q'|| | t/i ±0.1l| extracted in where the B — > ir form factors were computed 

in the light-cone sum rule (LCSR). The distinction can be traced back to the different q 2 dependence of the form 
factor f + {q 2 ) predicted in the kf factorization and in LCSR, albeit with the similar / + (0) value in both approaches. 
More dedicated efforts on the study of the shape of B — > ir form factors in QCD is in demand in order to resolve the 
potential difference in the extraction of \V u b\. 
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IV. CONCLUSION 

In this paper we have calculated the NLO corrections to the B — > n transition form factors at leading twist in 
the kr factorization theorem. Both the collinear and soft divergences in the NLO quark diagrams and in the NLO 
effective diagrams for meson wave functions are regularized by the off-shellness k\ of light partons. The b quark 
remains on-shell, such that it can be approximated by the standard effective heavy quark in the kr factorization. The 
key is that the soft gluons radiated by the b quark and attaching to other particle lines can be regularized by the 
virtuality of other particle lines. The NLO pion wave function is the same as constructed in the pion transition and 
electromagnetic form factors, confirming its universality. Compared to the pion wave function, the NLO B meson 
wave function contains the additional double logarithm \t? ((I / m 2 B ) . Because of the assumed hierarchy £i ^> m%, the 
appearance of this double logarithm demands the implementation of the resummation technique, which is expected 
to minimize the scheme dependence from different choices of £1. This subject, together with the asymptotic behavior 
of the B meson wave function in the kr factorization, will be discussed in a forthcoming work. 

The exact cancellation of the infrared divergences between the quark diagrams and the effective diagrams verifies 
the validity of the kr factorization for the B meson semileptonic decays at NLO level. Though the NLO hard kernel 
for the B — > n transition form factors contains a huge constant term, it is reduced by the large double logarithm 
hv ((^ / m 2 B ) mentioned above. This is the reason why the conventional choice of the factorization scale in the PQCD 
approach, as the virtuality of internal particles, can work to render the NLO corrections under control. By tuning 
the renormalization scale to cancel the single-logarithmic and constant terms, which is still lower than the B meson 
mass in the dominant kinematic region, the NLO corrections are about 30% of the form factors. The effect of varying 
the meson wave functions has been also investigated: the model for the B meson wave function with a stronger 
suppression at a small momentum fraction, and the asymptotic model for the pion wave function lower the NLO 
corrections down to 20%. 
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